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1 Introduction 



The models of relativistic particles with spin have been proposed for a long 
time. The first published work concerning the lagrangian description of the 
relativistic particle with spin was the paper by Frenkel which appeared in 
1926 [1]. After that the literature on the particle with spin grew vast [2]. 

The models involving only conventional coordinates are called the classi- 
cal models while the models involving anticommuting coordinates are gener- 
ally called pseudo-classical. 

In this paper we shall confine ourselves to discuss the relativistic spin one 
half particle models involving anticommuting vectorial degrees of freedom 
which are usually called the spinning particles. Spinning particles are in some 
sense the classical limit of the Dirac particles. After the first quantization 
these new anticommuting variables are mapped into the Dirac matrices and 
they disappear from the theory [3,4]. 

The action of spin one half relativistic particle with spinning degrees of 
freedom described by Grassmannian (odd) variables was first proposed by 
Berezin and Marinov [5] and soon after that was discussed and investigated 
in the papers [6-10]. 

In spite of the fact that the anticommuting Grassmann variables do not 
admit a direct classical interpretation, the lagrangians of these models turn 
out to be suitable for the path integral description of the quantum dynamics. 
The pseudo-classical equations acquire physical meaning when averaged over 
inside the functional integrals [5,11]. In the semi-classical regime, neglecting 
higher order quantum correlations, it should be allowed to replace some com- 
binations of Grassmann spin variables by real numbers. Using these ideas 
the motion of spinning particles in external fields have been studied in Refs. 
[5, 12-14]. 

On the other hand, generalizations of Riemannian geometry based on 
anticommuting variables have been proved to be of mathematical interest. 
Therefore the study of the motion of the spinning particles in curved space- 
time is well motivated. 

In the present paper we investigate the motion of pseudo-classical spin- 
ning point particles in curved spaces. The generalized Killing equations for 
the configuration space of spinning particles (spinning space) are analysed 
and the solutions of the homogeneous part of these equations are expressed 
in terms of Killing- Yano tensors. We mention that the existence of a Killing- 
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Yano tensor is both a necessary and a sufficient condition for the existence 
of a new supersymmetry for the spinning space [15-17]. 

The general results are applied to the case of the four- dimensional Eu- 
clidean Taub-NUT spinning space. The motivation to carry out this exam- 
ple is twofold. First of all, in the Taub-NUT geometry there are known 
to exist four Killing- Yano tensors [18]. From this point of view the spin- 
ning Taub-NUT space is an exceedingly interesting space to exemplify the 
effective construction of all conserved quantities in terms of geometric ones, 
namely Killing- Yano tensors. On the other hand, the Taub-NUT geometry 
is involved in many modern studies in physics. For example the Kaluza- 
Klein monopole of Gross and Perry [19] and of Sorkin [20] was obtained 
by embedding the Taub-NUT gravitational instanton into five-dimensional 
Kaluza-Klein theory. Remarkably the same object has re-emerged in the 
study of monopole scattering. In the long distance limit, neglecting radia- 
tion, the relative motion of slow Bogomolny-Prasad-Sommerfield monopoles 
is described by the geodesies of this space [21,22]. The dynamics of well- 
separated monopoles is completely soluble and has a Kepler type symmetry 
[18,23-25]. The problem of geodesic motion in this metric has therefore its 
own interest, independently of monopole scattering. 

The plan of this paper is as follows. In Sec. 2 we summarize the rel- 
evant equations for the motions of spinning points in curved spaces. The 
generalized Killing equations for spinning spaces are investigated and the 
constants of motion are derived in terms of the solutions of these equations. 
In Sec. 3 we analyse the motion of pseudo-classical spinning particles in the 
Euclidean Taub-NUT space. We examine the generalized Killing equations 
for this spinning space and describe the derivation of the constants of motion 
in terms of the Killing- Yano tensors. In Sect. 4 we solve the equations given 
in the previous Section for the special case of motion on a cone. This case 
represents an extension of the scalar particle motions in the usual Taub-NUT 
space in which the orbits are conic sections [18,23-25]. An explicit exact so- 
lution is given and, in spite of its simplicity, this solution is far from trivial. 
Our comments and concluding remarks are presented in Sec. 5. 
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2 Spinning spaces and Killing equations 



Spinning particles, such as Dirac fermions, can be described by pseudo- 
classical mechanics models involving anticommuting c-numbers for the spin 
degrees of freedom. The configuration space of spinning particles (spinning 
space) is an extension of an ordinary Riemannian manifold, parametrized by 
local coordinates {x^ 1 }, to a graded manifold parametrized by local coordi- 
nates {x 11 , ip^ 1 }, with the first set of variables being Grassmann-even (com- 
muting) and the second set Grassmann-odd (anticommuting) [3-17]. 

The dynamics of spinning point-particles in a curved space-time is de- 
scribed by the one-dimensional a-model with the action: 

S = £ dT (l x»x»+ l - g, v {x) V ■ (1) 

Here and in the following, the overdot denotes an ordinary proper-time 
derivative d/dr, whilst the covariant derivative of ip^ is defined by 

^ = r+x"rir- (2) 

The trajectories, which make the action stationary under arbitrary vari- 
ations 8x^ and Sip^ vanishing at the end points, are given by: 

^ = F + rtx A / = if^U/, (3) 

The antisymmetric tensor 

S^ = -i^ v (5) 

can formally be regarded as the spin polarization tensor of the particle. The 
first equation of motion (3) implies the existence of a spin dependent gravi- 
tational force [14] 

-n^ = l SKXR ^ ( 6 ) 

which is analogous to the electromagnetic force, with spin replacing the elec- 
tric charge as the coupling constant. The second equation of motion (4) 
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can be expressed in terms of this tensor (5) and it asserts that the spin is 
covariantly constant 

The interpretation of S^ u as spin tensor is corroborated by studying elec- 
tromagnetic interaction of the particle [5,9,13,14]. From such an analysis 
it results that the space-like components are proportional to the magnetic 
dipole moment of the particle, whilst the time-like components S° l represent 
the electric dipole moment. The requirement that for free Dirac particles 
the electric dipole moment vanishes in the rest frame can be written as a 
covariant constraint [3] 

g uX S^x x = (8) 
which, in terms of the Grassmann coordinates, it is equivalent to 

= 0. (9) 

The concept of Killing vector can be generalized to the case of spinning 
manifolds. For this purpose we consider the world-line hamiltonian given by 

H = \g^^ v (10) 

where 

n M = g^x" (11) 

is the covariant momentum. 

For any constant of motion J(x, IT, tp), the bracket with H vanishes 

{H,J} = (12) 

where the Poisson-Dirac brackets for functions of the covariant phase space 
variables (x, H, ip) is defined by 

The notations used are 

BF BF i 
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and ap is the Grassmann parity of F : of = (0, 1) for F=(even,odd). 
If we expand J7"(x, II, ip) in a power series in the covariant momentum 

CO 1 

J = E -^-^(x, von M1 . . . n Mn (is) 

n=0 n - 

then the bracket {H, J} vanishes for arbitrary 11^ if and only if the compo- 
nents of J satisfy the generalized Killing equations [3,15,26] : 

d 7^ 

J (Hi...H n ;ii n +i) + ^ Sn+OA^ - 2 V V ^pcv^^J ^... Mn ) 

where the parentheses denote symmetrization with norm one over the indices 
enclosed. 

In the scalar case, neglecting the Grassmann variables {ip 1 *}, all the gen- 
eralized Killing equations (16) are homogeneous and decoupled. The first 
equation (n = 0) shows that is a trivial constant, the next one (n = 1) 
is the equation for the Killing vectors and so on. In general, for a given 
n, neglecting all spin degrees of freedom, eq.(16) defines a Killing tensor of 
valence n 

^..^i>(*) = ° ( 17 ) 

and from eq.(15) 

J = J^(x)U^ . . .U^ (18) 

is a first integral of the geodesic equation [27]. 

In the spinning case the symmetries can be divided into two classes. First, 
there are four independent generic symmetries which exist in any theory : 

1. Proper-time translations generated by the hamiltonian H (10); 

2. Supersymmetry generated by the supercharge 

Qo = IW; (19) 

3. Chiral symmetry generated by the chiral charge 

r, = ^^..^---^ w ; (2°) 
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4. Dual supersymmetry generated by the dual supercharge 

Q* = *{r*,Qo} = (rf^Vfo*...«n'»ir ■■■V d (21) 

where d is the dimension of space-time. 

As a rule we have the freedom to choose the value of the supercharge Q 
and any choice gives a consistent model. The condition for the absence of an 
intrinsic electric dipole moment of physical fermions (leptons and quarks) as 
formulated in eq.(9) implies 

Q = 0. (22) 

However, for the time being, we shall not fix the value of the supercharge, 
keeping the presentation as general as possible. 

The second kind of conserved quantities, called non-generic, depend on 
the explicit form of the metric g^x). In the recent literature there are exhib- 
ited the constants of motion in the Schwarzschild [28], Taub-NUT [17,29-33], 
Kerr-Newman [15,16] spinning spaces. 

In what follows we shall deal with the non-generic constants of motion in 
connection with the Killing eqs.(16) looking for the general features of the 
solutions. The spinning particle constants of motion can be seen either as 
extensions of the constants from the scalar case or new ones depending on 
the Grassmann- valued spin variables {ip^}- 

Let us assume that the number of terms in the series (15) is finite. That 
means that, for a given n, vanishes and the last non-trivial gener- 

alized Killing equation from the system (16) is in fact homogeneous : 

Qj{n) 

The line of action to solve the system of coupled differential equations (16) 
is standard. One starts with a solution of the homogeneous eq.(23) 

which has to be introduced in the right-hand side (RHS) of the generalized 
Killing equations (16) for jj^~^ n l and the iteration must be carried on to 
n = 0. 

For the beginning let us note that eq.(23) has solutions which do not 
depend on the Grassmann coordinates. These are the Killing tensors of 
valence n, as it can be seen comparing eq.(17) with eq.(23) in which all spin 
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degrees of freedom are neglected. However, for the spinning particles, the 
generalized Killing equations (16) are not decoupled. Even if one starts with 
a Killing tensor of valence n we get from the remaining Killing equations the 
components Jjj^\ m < n) with non-trivial spin contributions. 

Therefore the quantity (18) is no more conserved and the actual constant 
of motion is 

n i 

J=H -J^.L m m ■ ■ ■ (24) 

m=0 m - 

in which J'j^} -ftm (x, ip) with m < n has a non-trivial spin dependent expres- 
sion. 

The construction of the conserved quantity (24) in which the last term 
^ n is a Killing tensor can be done effectively. We shall illustrate this 
construction with a few examples. For n = eq.(17) is satisfied by a simple, 
irrelevant constant. The first non-trivial case is n — 1. In this case eq.(17) 
is satisfied by a Killing vector R^. 

R M = 0. (25) 

Introducing this Killing vector in the RHS of the generalized Killing 
eq.(16) for n = one obtains for the the expression [17]: 

J {0) = \R\,-aVV (26) 

where the square bracket denotes antisymmetrization with norm one. Con- 
sequently, starting with a Killing vector we get in the spinning case the 
conserved quantity (24) in the form 

J = \RwV , P + RJV t . (27) 

A more involved example is given by a Killing tensor jff = satisfy- 
ing eq.(17) for n = 2: 

= 0. (28) 

This solution must be introduced in the RHS of the generalized Killing 
equation (16) for and then we have to evaluate the new JM. Unfor- 
tunately it is not possible to find closed, analytic expressions for J"« and 
j7"(o) involving the spin variables using directly the components of the Killing 
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tensor K^ v . But assuming that the Killing tensor K^ u can be written as 
a symmetrized product of two Killing- Yano tensors, the construction of the 
conserved quantity (24) is feasible. 

We remind that a tensor is called a Killing- Yano tensor of valence 

r [27,34] if it is totally antisymmetric and it satisfies the equation 

/mi-^-i(av;A) = 0. (29) 

It is known that the Killing- Yano tensors play a key role in the Dirac 
theory on a curved space-time [35]. The study of the generalized Killing 
equations strengthens the connection of the Killing- Yano tensors with the 
super symmetric classical and quantum mechanics on curved manifolds. 

For the generality, let us assume that the Killing tensor K^ u can be written 
as a symmetrized product of two different Killing- Yano tensors 

Kg = \(ftxfj X + fixft) (3°) 

where is a Killing Yano tensor of valence 2 and type i. We use for the 
Killing tensor two additional indices i,j to emphasize the fact that it is 
formed from two different Killing- Yano tensors. 

Introducing the Killing tensor K^ v in the form (30) in the RHS of eq.(16) 
for n — 1 we can express the solution J^ 1 ^ in terms of the Killing- Yano 
tensors and their derivatives [15,32]: 

= y> x r(f?,D„f? x + fr a D v f? x 

+ \frc J x„ P + fYc lXap ) (31) 

where the tensor c ifiu \ is 

CitivX = —1fi[v\;ii\- (32) 

Finally, using j£ )fl in the RHS of eq.(16) for n = we get for J^f 

= -^VVV T (iW//Vir + \ciXaCjpr,). (33) 

Collecting the quantities (30), (31) and (33) in eq.(24) we get the corre- 
sponding conserved quantity: 

Jij = ^K^U,U V + J^U, + 4 0) . (34) 
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Higher orders of the generalized Killing eq.(16) can be treated similarly, 
but the corresponding expressions are quite involved. On the other hand, for 
practical purposes (see Section 3), it turns out to be sufficient to consider in 
detail the Killing tensors up to the valence 2. 

In what follows we shall return to the eq.(23) looking for solutions de- 
pending on the Grassmann variables {t/^}. The existence of such kind of 
solutions of the Killing equation is one of the specific features of the spinning 
particle models. 

Even the lowest order eq.(23) with n = has a non-trivial solution [17,31] 

j {0) = \f»vW (35) 

where is a Killing- Yano tensor covariantly constant. Moreover, from 
eq.(24), we infer that is a separately conserved quantity. 

The next eq.(23) with n = 1 can have different kinds of solutions. The 
most remarkable class of solutions is represented by 

j$ = (36) 

generated from a Killing- Yano tensor of valence r. Again, introducing this 
quantity in the RHS of eq.(16) for n = we get for J^: 

j {0) = ^(-ir + %,.,^ +1 ] -r 1 - (37) 

and the constant of motion corresponding to these solutions of the generalized 
Killing equations is [31]: 

Q f = f^n^r 2 ---^ 

+^(-i) r+ %,.,v^ +1 ] • r 1 ■ (38) 

This quantity is a super invariant 

{Qf,Qo} = (39) 

for the bracket defined by eq.(13). A similar result was obtained in ref.[36] 
in which it is discussed the role of the generalized Killing- Yano tensors, with 
the framework extended to include electromagnetic interactions. This result 



10 



extends the analysis from refs. [15-17] where it is established that the existence 
of a Killing- Yano tensor of the usual type (r = 2) is both a necessary and 
a sufficient condition for the existence of a new supersymmetry of the type 
(38) obeying the superinvariance condition (39). 

To conclude, we mention that eq.(23) with n = 1 has also many other 
solutions by forming combinations of different Killing- Yano tensors. As a 
rule, the corresponding constants of motion are not completely new and 
they can be expressed in terms of the quantities described above. We shall 
illustrate this fact choosing for the solution of eq.(23) with n — 1 the quantity 

J™ = R»M x r (40) 

where is a Killing vector (Killing- Yano tensor with n — 1) and f\ a is 
a Killing- Yano tensor covariantly constant. Introducing this solution in the 
RHS of the eq.(16) with n — 0, after some calculations, we get for 
[31,32]: 

J {0) = '-R^fxaVr^V- (41) 

Combining eqs.(40) and (41) with the aid of eq.(24) we get the constant 
of motion : 

J = M»V (#aIT a + l -R [x ,^ x r) ■ (42) 

As expected, we recognize in this expression the conserved quantities (27) 
and (35). 

In the next Section we shall apply these general results concerning the 
solutions of the Killing equations for spinning spaces to the case of the four- 
dimensional Euclidean Taub-NUT spinning space. 
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3 EUCLIDEAN TAUB-NUT SPINNING 
SPACE 

The Kaluza-Klein monopole [19,20] was obtained by embedding the Taub- 
NUT gravitational instanton into five-dimensional theory, adding the time 
coordinate in a trivial way. Its line element is expressed as: 

ds\ = —dt 2 + ds 2 

= -dt 2 + V- 1 (r) [dr 2 + r 2 d6 2 + r 2 sin 2 9 dip 2 } 

+V(r)[dx 5 + A(f)df} 2 (43) 

where r denotes a three- vector r = (r, 9, <p) and the gauge field A is that of 
a monopole 

A r = A d = 0, A v = 4m(l - cos6>) 

B = rot A = — r . (44) 



The function V(r) is 



1 • — ) ' (45) 



r / 



and the so called NUT singularity is absent if x 5 is periodic with period lGnm 
[37]. 

It is convenient to make the coordinate transformation 

4m(x + y?) = -x 5 (46) 

with < x < 4tt, which converts the four-dimensional line element ds± into 

ds\ = V-\r) [dr 2 + r 2 d9 2 + r 2 sin 2 9 dip 2 } + lQm 2 V(r) [d X + cos 9 dip} 2 . (47) 

Spaces with a metric of the form given above have an isometry group 
SU(2) x 17(1). The four Killing vectors are 

D A = RfXd li , A = 0,1, 2, 3, (48) 
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where 



D 



o 



d_ 



^ d n d cos cp d 

Di = — sin ip — — — cos (p cot 9 — I ; — — — , 
a0 a<^ sinfe 1 ax 

d . d sin <£> 9 

L> 2 = cost/?— -sm^cote/ — + -— , 
a6> a<yC sine* ax 

1,3 = h < 49 > 

Do which generates the U(l) of x translations, commutes with the other 
Killing vectors. In turn the remaining three vectors, corresponding to the 
invariance of the metric (47) under spatial rotations (A = 1,2,3), obey an 
SU{2) algebra with 

[D 1 ,D 2 ] = -D 3 ,etc... (50) 

In the purely bosonic case these invariances would correspond to conser- 
vation of the so called "relative electric charge" and the angular momentum 
[18,23-25] : 

q = 16m 2 V(r) (x + cos 9 ip), (51) 
— » 

-» T 

j = fx p + q-. (52) 
r 

where p = V^ 1 (r)f is the "mechanical momentum" which is only part of the 
momentum canonically conjugate to f. 

As observed in [18], the Taub-NUT geometry also possesses four Killing- 
Yano tensors of valence 2. The first three are rather special: they are covari- 
antly constant (with vanishing field strength) 

Atyi 

fi = 8m(dx + cos 9d(p) A dxi — 6^(1 H )dxj A dxk, 

D,f7x = 0, * = 1,2,3. (53) 
Moreover, they are mutually anticommuting and square the minus unity: 

fif S + fifi = -2Sij- (54) 

Thus they are complex structures realizing the quaternion algebra. In- 
deed, the Taub-NUT manifold defined by (47) is hyper-Kahler and, as a 
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consequence, the corresponding supersymmetric a-model has an N = 4 su- 
persymmetry. 

The fourth Killing- Yano tensor is 

fy = 8m(dx + cos Odip) A dr 

+4r(r + 2m) (1 + sin Odd A dip (55) 

and has only one non-vanishing component of the field strength 

f Yre ^ = 2(1 + ^-)r sm6. (56) 

In the Taub-NUT case there is a conserved vector analogous to the Runge- 
Lenz vector of the Kepler-type problem whose existence is rather surprising 
in view of the complexity of the equations of motion. This conserved vector 
is: 



K = \K, V WW = p x j+ (£- - AmE^j P - 
where the conserved energy E, from eq. (10), is 



(57) 



E 



\r^ v = \v-\r) 



f-V 

\4m/ 



1 4m + r ., 1 -o 1 , , . 

r H — {4m + r)ru H — (4m + r) r sm #92 

2 



+8m 2 



4m + r 



(cos 6*0 + xf- 



(58) 



The components K ifiU involved with the Runge-Lenz type vector (57) are 
Killing tensors and they can be expressed as symmetrized products of the 
Killing- Yano tensors fa (53) and f Y (55) as in eq. (30) [17,18,33]: 



Ki^ u = m [fy^\fi X v + fYu\fi X ^) + — (Ro^Riu + RovRin) 



(59) 



This equation corrects some old formulas from the literature [18]. 

Using these conservation laws one can determine the orbits. Eq.(52) 
implies that 

— * 

-f T 



r 



\j \ cos 9 = q 



(60) 
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which fixes the relative motion to lie on a cone whose vertex is at the origin 
and whose axis is j. Moreover, taking into account the existence of the 
Runge-Lenz vector (57), one finds that the trajectories lie simultaneously on 
the cone (60) and also in the plane perpendicular to 

n = qlC + [AmE - £-^J j. (61) 

Thus they are conic sections. 

Starting with these results from the bosonic sector of the Taub-NUT 
space one can proceed with the spin contribution to the conserved quantities 
(51), (52) and (57). 

First of all, corresponding to the generic symmetries described in the 
previous Section, there are four universal conserved charges. For the Taub- 
NUT spinning space these are : 

1. The energy (58); 

2. The supercharge (19): 

Q = rip r + (Am + r)r9ib° 

r 



+ 



(Am + r)r sin 9 (f + q cos 9 tp v + qip x \ (62) 



3. The chiral charge 

r* = Am(Am + r)r sin 9 t/j r ip ^ ip x ; (63) 

4. The dual supercharge 

Q* = Am(Am + r)r sin 9(f ip^ i\) x — 9 i\) r ^ i\) x 

+ i pip r ^^ x -xV^ 9 ^)- (64) 

From eq.(4) which shows that ip^ is covariantly constant, we find that the 
rate of change of the spins is: 

• 2m . r 2 + 2mr ■ e 

W = —r A x r W +— ; 9w 

r[Am + r) Am + r 
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( r 2 + 2mr 0n 32m 3 rcos 2 #\ 

+ — A sm ° + ^ \T~ <fV 

\ 4m + r (4m + r) A J 

32mVcos^ . 32mV . 

(4m + ry (4m + r)^ 

;a r + 2m . , fl • ,_, 8mr + r 2 „ „ . ,, 

r(4m + rj (4m + rj 2 



(4m + r) 



2 



• r + 2m 8m 2 + 8mr + r 2 cos • 

^ = 71 — T( r ^ +^^) 71 — "T2 ^i 9 ^ 

r(4m + r) (4m + r) z smt/ 

(4m + r) 2 sm^ 

^ X = 77~~ r~ + , ^ t A r^ x + xV) 

(4m + r) r(Am + r) 

/ 8m 2 + 8mr + r 2 cos 2 6> 1 . \ .• . 

+ 77—TT2 77 + .^ ^ + ^ 

\ (4m + r) 2 sine* 2 y 

8m 2 COS6 1 . • , v . , fl . . . 

(#^ X + X^). (65) 



(4m + r) 2 sin 6* 

As a rule, the complicated eqs. (3) and (4) should be integrated to obtain 
the full solution of the equations of motion for the usual coordinates 
and Grassmann coordinates {V^}- I n addition to the brute force method of 
trying to solve these equations there are the general prescriptions described 
in the previous Section which are considerable simpler. In what follows we 
shall use the Killing- Yano tensors to generate the constants of motion for 
spinning particles. 

We start with the observation that the angular momentum (52) and the 
"relative electric charge" (51) are constructed with the aid of the Killing 
vectors (49). The corresponding conserved quantities in the spinning case 
are the followings: 

J = B + l (66) 
J = B + q (67) 
where we used eq.(27), and we introduced the notation: J = (Ji, J 2 , J3), B = 
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(Bi, B 2 , B 3 ). From eq.(26), the scalars Ba 



B a = % -Ra[w\W 



have the following detailed expressions: 



B 



Bo 



B* = 



+ 
+ 



+ 
+ 



+ 



32m 3 cos 9 
(4m + r) 2 



S rip + 



32m 3 _^ rx _ 8m 2 r sinfl ^, 



(4m + rf 



4m + r 



f , n ,^0 8m 2 rsin6> 
- sin ip\ (2m + r)S r& H : S vx 



cos (/? 



32m 3 



(4m + r) 2 



4m + r 
(2m + r) sin 6* cos 9S rif 



32m 3 sin ^ rv 8m 2 r + (8mr 2 + r 3 ) sin 6> 2 

rO x H 



(4m + r) 2 
8m 2 r cos 9 
4m + r 



4m + r 



// . wr9 8m 2 rsinfl 
cos</? (2m + r)<S H : b vx 



+ sin if 



32m 3 



(4m + r) s 



32m 3 sin 9 
(4m + r) 2 
8m 2 r cos 
4m + r 



4m + r 

(2m + r)^J sin cos 9S r(p 

8m 2 r + (8mr 2 + r 3 ) sin# 2 
4m + r 



(2m + r) sin6» 2 + 



32m 3 cos# 2 



(4m + r) 2 



^ 32m 3 cosfl ^ 



3mr 



r 3 ) sin cos 



3m 2 r 



(4m + r) 2 



4m + r 4m + r 

We mention that the above constants of motion are superinvariant: 

{Ja,Qo} = , A = 0,...,3. 
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Also, the components (66) of the angular momentum satisfy, as expected, 
the SO (3) algebra: 

{Ji, Jj} = e ijk J k , i,j,k = 1,2, 3. (71) 

We consider now the Killing- Yano tensors of valence 2 and we search for 
those constants of motion built of them. 

Using eq.(38) we can construct from the Killing- Yano tensors (53) and 
(55) the supercharges Qi and Qy. The supercharges Qi together Q from 
eq.(19) realize the N = 4 supersymmetry algebra [17]: 

{Qa,Qb} = -2i5 ab H , A,B = 0,...,3 (72) 

making manifest the link between the existence of the Killing- Yano tensors 
(53) and the hyper- Kahler geometry of the Taub-NUT manifold. Moreover, 
the supercharges Qi transform as vectors at spatial rotations 

{Qi,Jj} = €i jk Qk , i,3,k = 1,2,3 (73) 

while Qy and Qo behave as scalars. 

We note also that the bracket of Qy with itself can be expressed in terms 
of the hamiltonian, angular momentum and "relative electric charge" : 

{Qy, Qy} = -2, (H + • ( 74 ) 

On the other hand, the existence of the Killing- Yano covariantly constant 
tensors f\ (53) is connected with three new constants of motion as shown in 
eq.(35): 

Si= l -U„rr , i = 1,2,3 (75) 

which realize an 5*0(3) Lie-algebra similar to that of the angular momentum 
(71): 

{Si,Sj} = e ijk S k , %, j,k = 1,2,3. (76) 

These components of the spin are separately conserved and can be com- 
bined with the angular momentum J to define a new improved form of the 
angular momentum Jj = Jj — Si with the property that it preserves the 
algebra [17]: 

{Ii, Ij} = e ijk I k , i,j, k = 1,2,3 (77) 
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and that it commutes with the SO (3) algebra generated by the spin Si 



{I»S j } = 0. (78) 
Let us note also the following Dirac brackets of Si with supercharges 

{Si,Qo} = — 

{Si,Qj} = ^(SijQo + e ijk Q k ). (79) 

To get the spin correction to the Runge-Lenz vector (57) it is necessary to 
investigate the generalized Killing eqs.(16) for n — 1 with the Killing tensor 

— * 

in the RHS. For an analytic expression of the solution of this equation 
we shall use the decomposition (59) of the Killing tensor K^ v in terms of 
Killing- Yano tensors. Starting with this decomposition of the Runge-Lenz 
vector K from the scalar case, it is possible to express the corresponding 
conserved quantity K, in the spinning case [33]: 

JCi = 2m (-i{Q Y , Qi} + ^Wo) (80) 

This expression differs from previous results presented in the literature 
[17,31] and the difference has the origin in the corrected form of relation (59). 
A detailed expression of the components K.^ is: 



/G = 2m 



+ (^fi X /3fYfia;X + f^^Ya^X 



(RiatfRon + RoatfRiii) ) S al3 W 



16m 2 



More explicitly, we can write the Runge-Lenz vector K, for the spinning 
case as in eq.(34): 

K = ^IPIT + S^W + S (82) 
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where the first term is the Runge-Lenz vector (57) from the scalar case and 
the last two terms represent the specific spin contribution. A detailed ex- 
pression of this contribution is [30,33]: 



Si, ■ IP 



— (4m + r) cos 9 cos cp ■ S r0 
+ Amr cos 9 sin ip ■ 

S 6<P + 4mr 

sin ip ■ 

+ (4m + r) sin 9 sin ip ■ S rLp + 4mr sin 9 cos 9 cos ip • S^*] f 



+ 



tA \ ■ n nre 4mr(6m + r) . 
r(4m + r) sin 9 cos ip • S ; cos 9 sin ip • S * 



4m + r 



4mr(6m + r) 



sirup ■ S rx 



4m + r 

+r 2 (6m + r) sin 9 sin (^S* ^ — 4mr 2 sin 2 9 cos <^S"^ X 



+ 



4mr(6m + r) fl 

cos 9 sin 

4m + r 



+ r(4m + r) sin 3 # cos <p + — — sin 9 cos 2 cos ip S Tlf 



— Amr + 



8m 2 r 



4m + r (4m + r 



(4m + r) 3 

sin # cos 9 cos ^S"** 



+ 



+ 



V 2 (32m 3 + 64m 2 r + 14mr 2 + r 3 ) . 
(4m + r) 2 



sin 2 9 cos cos (p 



3„,2 



32m r 
(4m + r) 2 



cos 3 cos ^ 



+ ( 4mr 2 sin 2 cos <z> + r cog 2 ^ cog _ j g6»x 

\ (4m + r) 2 J 



3^,2 



32m t 
(4m + r) 



■ sin cos 9 sin <^S" PX 



+ 



4mr(6m + r) . e 256m 4 r . n rirv 
sm</?5 + T-; sm e* cos </?5 x 



4m + r 
+ I 4mr + 



m 2 r 



+ 



4m + r (4m + r) 3 



(4m + r) 3 

sin 9 cos cos ^S""^ 
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$1 



+ 



32m 3 r 2 



cos 2 cos ip 



v (4m + r) 2 

64m 3 r 2 + 16m 2 r 3 + 2mr 4 . 



(4m + r) 2 



sin 2 9 cos ^ 



32m 3 r 2 „ , 32m 3 r 2 . „ . 
+ — : — cos 9 cos ipb x — — — sin sin ^ x 



(4m + r)< 



(4m + r) s 



8m 2 r 2 (8m + r) sin 9 2 cos 99 e 

Ta \~5 ^ ^ ' 

(4m + ry 



X 



(83) 



(4m + r) sin 0S re - 4mr sin 2 9S* X 



+ 
+ 

+ 



r{Am + r) cos0S r9 - 4mr 2 sin 9 cos 9S m 



256m 4 r 
(4m + r) 3 
256m 4 r 



cos 3 + r(4m + r) sin 2 cos 9) S rtf 



cos 2 + 



4mr(6m + r) 



(4m + r) 3 
— r 2 (6m + r) sin 3 + 



4m + r 
g 6m 3 r 2 

(4m + r) 2 



sin 2 9 S rx 



sin cos 2 



+ 2mr 2 2 



16m 2 



(4m + r) 2 / 



sinflcosflS** 



95 



+ 



— 4mr sin — 



8m 2 r 
4m + r 



sin 2 + 



256m 4 r 
(4m + r) 



cos 



256m 4 r 32m 3 r 2 . 
+ —. — cos9S x — — — sm 9S x 



(Am + r) 3 (4m + r) 2 

24m 2 



- 4mr 1 + 



(4m + r) 2 



sin 9 cos 0S^ 



8m 2 r 2 (8m + r) cos sin e 



(84) 

(4m + r) 3 

The components S^IP and S 2 can be obtained from Si^rP respectively 
Si with the substitutions: 



sm</? 
cos ip 



COS if, 

sin </? . 



(85) 
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Therefore, in the spinning case, the Runge-Lenz vector contains addi- 
tional terms linear and quadratic in the spin. The presence of a contribution 
quadratic in the spin, non-existent in Refs. [17,30], is again related to the 
term JjJ from eq.(80). We would like to emphasize that the contribution 
of the JjJo term in eq.(80) is essential in reproducing the known vectorial 
expression of the Runge-Lenz vector in the scalar case, and gives the correct 
Poisson-Dirac bracket between two components of the Runge-Lenz vector. 

The Dirac brackets involving the Runge-Lenz vector (80) are (after some 
algebra) : 



{/Q,Q } = 0, 



T 2 

Jo 2H 



16m 2 



(86) 



and they are similar to those known from the scalar case. The Runge-Lenz 
vector JC together with the total angular momentum J generates an SO (4) or 
5*0(3, 1) algebra depending upon the sign of the quantity (3^2 — 2i?) | ^ M =o 
is positive or negative. 

In conclusion, all conserved quantities for motions in spinning Taub-NUT 
space have been expressed in terms of the geometric objects (49), (53) and 
(55). Other expressions involving the Killing- Yano tensors will produce con- 
served quantities which are not new, but rather combinations of the above 
primary conserved quantities. For example, let us consider a solution of the 
homogeneous eq.(23) for n = 1 of the type (40): 

j { 21 = iivwv, 

A = 0,...,3 , j = l,2,3. (87) 
After some algebra we get the constants of motion of the form (42): 

-AiSjJ A , A = 0,...,3 , J = 1,2,3. (88) 

As expected, the constants J A j are not new, being expressed in terms of 
the constants J a (62), (63) and Sj (75). However, the combinations (88) arise 



22 



in a natural way as solutions of the generalized Killing equations and appear 
only in the spinning case. Moreover, we can form a sort of Runge-Lenz vector 
involving only Grassmann components: 

Li = —tijkSjJk , k = 1,2,3 (89) 
m 

with the commutation relations like in eqs.(86): 

{Li, Jj\ GijfcLk, 

{L t ,Lj\ = (sf-S 2 )-^e tjk J k . (90) 
Note also the following Dirac brackets of Lj with supercharges: 

{Li,Q } = —-^tijkQjJk, (91) 



{Li,Qj} = ^(e ijk QoJk - SijQkMj; + QiMj). (92) 
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4 SPECIAL SOLUTION 



In spite of the fact that all conserved quantities have been expressed in terms 
of geometric ones in a close form, their detailed expressions (69), (83), (84) 
are quite intricate. 

We wish to consider a special class of solution of the equations of motion 
which is very simple, but not at all trivial. In this Section we confine our- 
selves to the motion on a cone on the analogy of the scalar case where the 
trajectories are conic sections. 

For this purpose let us choose the z axis along J so that the motion of 
the particle may be conveniently described in terms of polar coordinates 

r = re(9,<p) (93) 

with 

e = (sin 6 cos (p, sin 9 sin ip, cos 9) . (94) 
For this choice of the axis we have: 

(4m + r)r9 = -(2m + r)S re - J^!ll sin QS m , (95) 

4m + r 



<P = 



+ 
+ 



g 

r(4m + r) cos# 
32m 3 



32m 3 



2m + r 



r(Am + r) 3 r(4m + r) 



r(4m + r) 3 cos# 
8m 2 



S r X + 



8m 2 



f8mr + r 2 ) sin# 2 



(4m + r) 2 sin 9 cos 9 



grip 



s e<p 



{Am + r) 2 sin 9 
and from eqs.(66), (67) and (69) 

Jf_ 
J — 



s e x 



- r (4m + r)S ev sin 9. 



(96) 



(97) 



In what follows we shall consider the angle 9 = constant, i.e. 9 = 0. The 
rate of change of spin (65) can be written in terms of the antisymmetric spin 
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polarization tensor S^ v (5) as follows: 



S r0 = - S re + J™ + \, qsme S r * 

4m + r 2(4m + r) 3 

8m 2 q sin 9 grx r sin 2 9 + 2m q g9ip 



r(4m + r) 3 cos# (4m + r) 2 cos# 

2mg 



(4m + r) 2 



ory 5 orS ^ rry 2mq Q^PX 

2r{4m + r) sin 4m + r (4m + r)" 1 

grx _ ? Cjr8 _|_ ^ COS ^ gr<£ 

2r(4m + r) sin cos # 4m + r 

(r sin 2 9 + 2m)q ^ 
(4m + r) 2 COS6 1 

^ = (2m + r)g grf) g 2m + r ^ ^ 
r 2 (4m + r) 2 cos # r(4m + r) 

8m 2 sinflg gyx 
r(4m + r) 3 cos 6* 

^ = — — s r9 + —!—coses e v--s e x 

8m(4m + r) 2 4m + r r 

8m + r • a era 
2(4m + r) 3 

gipx 1 g r f (2m + r) q grx 

8m(4m + r) 2 r 2 (4m + r) 2 cos # 

9 ^ Q 



$o<p 1 s&x 



2r(4m + r) sin # cos # 2r(4m + r)sin# 



r 
r 



Since we are looking for solutions with 9 = we have from eq.(95) 

(2m + r)(4m + r) 
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This relation implies that the special solution investigated in this Section 
is situated in the sector with 

r* = o. (ioo) 

We mention that in this sector the system of eqs.(98) is satisfied even if 
the angle 9 is not constant. 

Using eq.(99) we can express S r6 through S vx and the following equations 
are equivalent to the system (98) 



^-\(Am + r)S^} 
at 



— \cos9S r,p + S rx } 
at 



d r 
di 



r{Am + r)S 9lp ] = -2 



Am + r 
0, 

sin 9 



qS m , 



d_ 
di 



r cos 9 S e * + rS dx 



cos 9 Am + r 
sin 9 r 



qS**, 



4m 4m + r 



qS vx . 



(101) 



Thus the equations of motion for S^ v are written in a more tractable form 
and the solution follows without difficulties. If we take into consideration the 
constraint coming from eq.(97), namely: 



d_ 
~dt 



r(Am + r)S e * 







then we have to impose 

q ■ S* x = 0. 

on the system (101). We shall analyze both solutions S vx 
successively. 

For S vx = 0, from eq.(99) we have also 



(102) 

(103) 
and q = 



S re = 0. (104) 
In spite of this drastic simplification, eqs.(101) have a non-trivial solution: 

(sin0=Fl) _ 



cos #(4m + r) 
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S r X - _ ( sin °T 1) E 

4m 4m + r 

S 0V = —r^ rE, 

r(4m + r) 

^ = -(d°fl±l)taiifl s _ cos^ 

4mr r(4m + r) 

where E is a Grassmann constant, commuting with if} 11 , and anticommuting 
with itself. 

In the case of this particular solution, from eqs.(67)-(69) we get that the 
spin contribution to the "relative electric charge" vanishes (B = 0) and 

Jo = q- (106) 

Therefore the "relative electric charge" has the same expression as in the 
scalar case. However the total angular momentum is modified by the spin 
contribution: 

j Q - — = q - Jcostf = -Esintf. (107) 

r 

Here J is the magnitude of the total angular momentum and eq.(107) fixes 
the angle 9 in terms of the constants q, J and E. Also the equations for ip 
and x are modified: 



<P 



r(4m + r)cos# (4m + r) 2 cos6*' 



8m + r E 
X = 16m 2 (4m + r) qT (4m + r) 2 ' ( ' 

At last, f can be derived from the energy, eq.(58). 

Concerning the second possibility, namely q = in eq.(103), we have from 
eqs.(lOl): 



Qrip 



4m + r ' 

cos 

4m + r 



r{Am + r) ' 
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r>6 x rs<f 
S e * = — - costf 



S 



r6 



r r(4m + r) ' 



4m + r ' 

S ra = , (109) 

r 

where C^ v are Grassmann constants of the same kind as S. 
Taking into account the constraints: 

8m 2 (4mC rx -sin 6C e f) ,,, n . 

J = B = ■ — , (110) 

(4m + r) z 



and, from (99), where we substitute (109): 

i 2 sin 9C 
2m + r 



c „ + 8"^m^ =0 (111) 



we get 

C r8 = Q ^ 0^ = (112) 

and 

C rx = , = (113) 

or 

C r X = S _^l C V (n4) 

4m v ' 

In conclusion, the case q = is included into the previous case (S^* = 0). 
Practically we must impose in addition to the previous solution the condition 
that q = and we get either the solution (112), (114) or a trivial one without 
any spin contribution (112), (113). 

Concerning the Runge-Lenz vector for 9 = constant we have from eqs.(83) 
and (84): 

S^IP = =f2 tan 9 sin ipfT, =p 2r tan 9 sin 2 9 cos (pipY, 
(8m + r) sin # cos (j? _ 
4m(4m + r) 
S 3/1 IP = ^2rsin 2 ^S 

tan #((8m + r) sin =F (4m + r)) _ . 

± "T^ 115 

4m(4m + r) 
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Again the component S^IP can be obtained from Si^IP with the substitu- 
tion (85). 

Using then: 

Sp= - (4m + r) " sin ^, (116) 
Avar 

pr (4m + r)r ni7^ 
~ ~y ' \ ') 

Jp= (g + £sinfl) — . (118) 

r 

we get to the conclusion that, for the case of motion lying on a cone, there 
is a conserved vector n orthogonal on p: 

-> (q — X sin 9) J ( 4 ^ q(q — X sin 9) \ /,,„s 

n = qK + '— AmE - ™— '- 119 

q \ Am J 

which is similar to the scalar case, eq.(61). Physical observables are obtained 
by averaging with some suitable density over the anticommuting parameters. 
After the integration was performed we may treat X as a classical variable. 
Therefore the trajectories of a spinning particle constrained to the motion 
on a cone are conic sections determined by the condition tip = 0. 
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5 CONCLUDING REMARKS 



The spinning particle model is a world line supersymmetric extension of the 
theory of a scalar particle. It describes a relativistic particle with spin one 
half. 

It is a theory which describes in a pseudo-classical way a Dirac particle 
moving in an arbitrary <i-dimensional space-time. In addition to the usual 
space-time coordinates, the model involves anticommuting vectorial coordi- 
nates which take into account the spin degrees of freedom. It is worth em- 
phasizing that along the world line of the particle there is a supersymmetry 
between the fermionic spin variables and the bosonic position coordinates. 
The model is a one-dimensional supersymmetric field theory on the world 
line. 

The term pseudo-classical refers to the fact that there is no classical 
interpretation for the anticommuting variables. To get the "observable" tra- 
jectories one has to average over the spin variables. On the other hand it is 
possible to quantize the model giving rise to supersymmetric quantum me- 
chanics. After quantization the conservation law for the supercharge becomes 
the Dirac equation. 

The constants of motion of a scalar particle in a curved space-time are 
determined by the symmetries of the manifold, i.e. if a space-time admits 
a Killing tensor K^...^ of valence r, then the quantity K jJil ... jJir U!l ■ ■ ■ IT^ is 
conserved along the geodesic. 

In the spinning case the generalized Killing equations (16) are more in- 
volved and new procedures should be conceived. The aim of this paper was 
to point out the important role of the Killing- Yano tensors to generate solu- 
tions of the generalized Killing equations. We presented a detailed discussion 
on how to construct conserved quantities out of Killing- Yano tensors for the 
Taub-NUT spinning space. Finally we solved the equation of motion for the 
case when the angle 9 is held fixed. This solution is most simple but far from 
trivial and the trajectories are the analogous of the ones of a scalar particle, 
being conic sections. 

The extension of these results for the motion of spinning particles in 
spaces with torsion [38] and/or interacting with background fields [3,36] will 
be discussed elsewhere [39]. 
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